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A spherically symmetric wormhole in Newtonian gravitation in curved space, enhanced 
with a connection between the mass density and the Ricci scalar, is presented. The 
wormhole, consisting of two connected asymptotically fiat regions, inhabits a spherically 
symmetric curved space. The gravitational potential, gravitational field and the pressure 
that supports the fluid that permeates the Newtonian wormhole are computed. Parti¬ 
cle dynamics and tidal effects in this geometry are studied. The possibility of having 
Newtonian black holes in this theory is sketched. 


1. Introduction 

The idea that Newton’s theory of gravitation, can be formulated in curved space 
has been recently analyzed by Abramowicz and collaborator£10. In addition to 
considering Newtonian gravitation in curved space, an equation linking the geom¬ 
etry of the 3-space with the matter sector was put forwardSEl. This is called an 
enhanced Newtonian gravitation. In particular, this modification to Newton’s the¬ 
ory allows the construction of a wormhole spac^. Wormholes in general relativity 
have been studied in several works se^^^®^for example. In this paper, we report on 
the construction of a spherically symmetric wormhol^ in this enhanced Newtonian 
gravitation. We further analyze the tidal effects that emerge from the gravitational 
forces and the curvature of space in this spherically symmetric Newtonian wormhole 
space. 

The outline of the paper is as follows: in Sec. we start by writing the equations 
that define the enhanced Newtonian theory of gravitation. In Sec. we construct a 
static spherically symmetric Newtonian wormhole, we find the gravitational field, 
the gravitational potential and the pressure of the fluid that supports the wormhole. 
Then, we study test particle’s motion in the wormhole geometry and gravitational 
field. Finally we analyze the tidal forces exerted by the gravitational field and the 
curvature of space in two nearby particles. In Sec.|^ we conclude and speculate on 
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the possible existence of truly Newtonian black holes in this enhanced Newtonian 
theory of gravitation. 

2. The fundamental equations 

The classical formulation of Newton’s theory of gravitation was constructed for an 
absolute 3-dimensional Euclidean space. We see, however, that the set of equations 
that comprise Newton’s theory of gravitation are also well defined for curved space. 
Indeed, Abramowicz et alEB recently proposed a formulation of Newton’s theory 
of gravitation in curved space. 

Poisson’s equation in curved space is given by 

= ( 1 ) 

where gij is the curved space metric, is the covariant derivative induced by the 
metric, the indices i,j run as i,j = 1,2,3, and p is the density of the matter. For 
static systems the continuity equation is trivially verified and the Euler equation is 
simplified to 

Vjp-kpVi0 = O, (2) 

where p is the fluid’s pressure. Now, a possible enhancement to Newton’s theory of 
gravitation was proposed irP where a relation between the geometry of space and 
the matter was introduced and given by 

R=2kp, (3) 

where R is the Ricci scalar and k is an arbitrary constant. The equations of motion 
of a test particle with mass m subjected to a gravitational potential </> are given by 
Newton’s second law, 

ma^ =—mg^^Vjfj), (4) 

where it was assumed that the inertial and gravitational mass of the test particle 
are equal. Eqs. 0-0 define an enhanced Newtonian gravitation. 

3. Spherical Newtonian wormholes 

3.1. Construction of the spherical Newtonian wormhole 

3.1.1. Matter density and metric 

Our purpose is to make use of Eqs. 0-0 to construct a static, spherical symmetric 
wormhole space. Let us then generically write the metric of the space as 

= A (r) dr^ + {d0^ + sin^ 9 dp^) , (5) 

such that the Ricci scalar R is 

2 [(A (r) - 1) A (r) -f rA' (r)] 


r^A (r) ^ 


( 6 ) 
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Now, to proceed we consider the expression for the mass density to be 

p{r)=ae-^ (7) 

where a has dimensions of mass density and b of distance. Substituting Eqs. and 
Q in Eq. ^ we find the general form of the metric, for the imposed mass density, 

(dd^ + sin^ 0 , (8) 

r + /3 re~^ + Ci 

where Ci is an integration constant and 


fj = b'^ka. 


(9) 


3.1.2. The embedding: the value of the integration constant and the 
parameter jd 


Having found the generic form of the metric, we have now to restrict the values of 
the parameters Ci , a and k to have a wormhole geometry. For such, we will follovEl 
and use an embedding diagram. So, in the Euclidean embedding space the axially 
symmetric embedded surface can be, using cylindrical coordinates (r, z, (p)^ uniquely 
described by a function z{r). Identifying the coordinates {r,ip) of the embedding 
Euclidean space with the coordinates (r, p) of the wormhole space, we find the 
following relation 


dz 

dr 


= ±. 


■ /3 re 


Cl 


- 1 , 


( 10 ) 


which can be used to study the properties of the space. From Eq. (10) and imposing 
the throat condition we find that the integration constant Ci is given by 


Cl = -b 




( 11 ) 


and the range of values of the parameter fd = kab"^ is constrained to be — oo < /3 < e. 
The same restriction to the parameter jd is found from the flare-out condition. 
Substituting Eq. (Ill in Eq. ([I]) we find 

1 


ds = 


1 _ Mr) 


-I- (dd^ -I- sin^ Odip'^) 


( 12 ) 


where 

b (r) = b + — ^b — re^~^^ . (13) 

Another restriction to have a wormhole geometry is that although near the worm¬ 
hole’s throat, r = b, the radial coordinate r is ill behaved, the proper radial distance 
must be finite. This condition implies a new restriction to fd 


^ ^ jd ^ jdcrit j 


( 14 ) 
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- — 1 

where the critical value of /3 is given by /3crit = infi< ^ < '' 

represents the infimum of the function in the specified range. This gives 

Pcrit = 2.338, 


, and inf 


(15) 


up to the third decimal place. Comparing Eq. (15) with the restriction found from 
the throat condition, /3 < e, we conclude that the product [} = kab^ is then re¬ 
stricted by Eqs. (141 and (15l. 

Due to the range of values that the parameter {3 might take, various distinct 
cases could be considered. Here we shall be interested in the class of wormholes 
whose parameters a, k and b obey 0 < ,5 ^ /3crit together with a > 0, and so fc > 0. 


3.1.3. Removal of the coordinate singularity 

Now that we have found the metric of the Newtonian wormhole we have to remove 
the coordinate singularity at r = 6. Considering a new coordinate I, defined as 
P = r'^ — 6^, such that —oo < I < oo, one finds that the metric of the space can be 
rewritten as 

ds^ = — - y - ^ ^-I- -I- 6^) -I- sin^ 0 dip'^) . 

{b'^ + P) \ e + /3e~^j — by/b'^ + P {e + P) 

(16) 

It will also be useful to write the mass density in the new coordinates. It is given 
by the expression 


p(0 


= a 


e~'^~^ (62 + 2P) 
b^ + P 


(17) 


3.2. Gravitational field, gravitational potential and pressure 
support of the spherical Newtonian wormhole 

Having defined the geometry of the space we can now study the gravitational po¬ 
tential of the Newtonian wormhole by solving Eq. Q. It is, however, considering 
the symmetries of the system, straightforward to solve Eq. Q by using Gauss’s law 
for gravity and integrating this equation over a volume whose boundary are the 
surfaces of constant gravitational potential. Defining the gravitational force field as 

= ( 18 ) 

with a comma denoting a partial derivative, we find 


(0 = 


Gm{l) 
2{b^ + P) 


\ 


(b^ + P) (e + Pe-^'^ -bVWTP{e + P) 


eP 


(19) 
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where the mass within radius I is given by 

pi-\-i^/b^ -X /2r — 11 

m{l) = 4:7Tab^ . dx. (20) 

^ Wx (1 +/3e-^) - (l + f) a/x 

Defining the magnitude of the gravitational field as 

g = ( 21 ) 

we give in Fig. its behavior as a function of the radial coordinate I for various 
values of the parameter p. 



Fig. 1. Magnitude of the gravitational force field for various values of the parameter /? = b^ka as 
a function of the radial coordinate I (here assume I > 0). When 0 = dcrit = 2.338 the magnitude 
of the gravitational force field tends to infinity at some value of l/b. 


Using the expression for the gravitational field, Eq. (191 together with Eq. (181, 
the gravitational potential (j) can be computed at a given point using numerical 
methods, provided the wormhole parameters are given. Moreover, these equations 
can also be used to find the pressure of the fluid that covers the wormhole space. 
Numerically solving Eqs. ([^ and (191 we find that the pressure is positive through¬ 
out the space showing that the wormhole is hold against gravitational collapse by 
pressur^. From Eq. ( [T^ one sees that the gradient of the potential (f) is the grav¬ 
itational force field 5®, and from Eq. ([^ the gradient of the pressure p goes with 
minus the gradient of the gravitational potential </>. Thus one can infer from Fig.j^ 
and the help of Eqs. (§. ( [T^ , and ( [M] ), that at the wormhole’s center I = 0, i.e., 
at the throat, the pressure is maximal, as one would expect. 

In the limiting case that the wormhole has P = /3crit a Newtonian event horizon 
develops at some radial coordinate Ih > 0, in the sense that any particle inside 
the sphere defined by > 0 can only affect an outside observer if it has infinite 
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acceleration, and thus infinite velocity. Thus, since no particles can come out of 
this inner region one is in the presence of a Newtonian black hole. However, for 
/3 = /3crit the pressure support at the throat goes to infinity, and so, in this limit, 
the wormhole eventually collapses. 


3.3. Test particle dynamics and motion in the spherical wormhole 
gravitational field 

3.3.1. The equations of motion 

Now we use Eq. Q to study the motion of a test particle in the gravitational field 
of the wormhole space. Consider that a particle’s path is described by a curve 7 , 
whose components are given by x® = {I (t) ,0(t),(p (t)) such that the particle’s 
velocity is u® = x®, where a dot means a derivative with respect to time t. The right 
hand side of Eq. Q is the gravitational field given in Eq. (19) and the left hand 
side the acceleration of the particle, defined as a® = Vj u®. It is possible to show 
that one can treat the problem of the motion of a test particle by considering pure 
equatorial orbits, i.e., 9 = ^,9 = 0 and 9 = 0. Gathering these results the equations 
of motion are 


1 + 


6(26^ + Z 2 ) (e + /3) - 2\/62 + P 

PP {l + P/b‘^)+b^ ( 

+P^ 

2Z(62 + Z 2 ) 

6 e^ {e + P) — y/P + Z 2 

+P 

)] 



+ 


by/Wfn? (e + /3) — ( 6 ^ + + /3 e 


d 


V+ 


= 


Gm{l) 


2sfe (52 + Z2) I 


+ (e 


9 = 9 = 0. 

2 


+ 2 


/ 


62 + ;2 


lip = 0. 


’e ^ ^ — 6\/52 + P {e + p), 

( 22 ) 

(23) 

(24) 


3.3.2. Solutions of the equations of motion 

Having found the equations of motion of a test particle, we have now to solve them. 
Let us consider the simpler problem of a test particle describing pure circular motion 
in the Newtonian wormhole space. In this case the coordinate I is a constant, Iq, 


say, so l = 0 and I = 0. From this, Eq. (22) can be solved for p and obtain 


ip = 


\ 


i/eGm (Zo) 


2(62 + Z2) J(62 + Z2) e + /3e--^ _5y6^(e + /3) 


(25) 
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Eq. (25) relates the radial position Iq and the angular velocity (p for a particle 


to describe a circular orbit in the wormhole space. Now, except the case of pure 


circular motion, to solve the equations of motion (22)-(241 for a more generic motion. 


numerical methods must be used, provided the initial position and velocity of the 
test particle and the wormhole parameters b and a are given. 


3.4. Tidal effects from the gravitational field and the geometry of 
the spherical Newtonian wormhole 

3.4.1. Tidal deviation equation for the spherical Newtonian wormhole 

We now analyze the relative separation of two test particles in the gravitational 
field of the spherical Newtonian wormhole space. The two test particles are initially 
considered to be infinitesimally close, subjected only to the gravitational field of 
the Newtonian wormhole, such that the acceleration of a test particle at a certain 
point is given by Eq. El and Eqs. (19l-(20). This relative separation, i.e., the tidal 


deviation equation, is in part given by the variation of the gravitational force field 
acted on each particle and another part given by the geodesic deviation due to the 
spatial curvature through the equation. 




= n^ i-NkG" 




(26) 


where n® is the separation vector between two infinitesimally close particles, is 


the gravitational force field generated by the wormhole mass, given by Eq. (19), u* 


is the velocity of the fiducial test particle and represents the Riemann tensor. 
Let us simplify the calculations, by remarking that, as was mentioned in 


Sec. (3.3), the velocity of the fiducial test particle can be assumed to be along 


the equator, 6 — it 12 (see Eq. (23)). Hence, in Eq. (26), we might take all the terms 
in to zero. Notice, however, that this does not mean that, in general, the accel¬ 
eration of the component n® of the separation vector is zero since the motion of the 
second particle may not be along the equator. Now, considering the metric of the 


spherical Newtonian wormhole space, Eq. (16), and assuming 0 = 7r/2 and 9 = 0, 


we find from Eq. (26) that the deviation equations for two infinitesimally close test 


particles in the spherically symmetric Newtonian wormhole space are 
D‘^n‘ 


dt'^ 

dt'^ 


= -n' {diG^ + T\iG^) + {n'^v^ - 




Riw W 


' R^ipo 





(27) 

{vy 

1 

(28) 



(29) 


where the Christoffel symbols that appear in Eqs. 


have the form 


^ u — 


+ e)ei^ (262 + f) - 2^6^ -h l^ 

Pp (62 -h P) + 6^ 


2bfi{b^ + P) 

b{j3 + e)e^ — Vb'^ + R 




(30) 



































8 Paulo Luz, Jose P. S. Lemos 


T^O _ ^ 

SI Vl b^ + p- 


(31) 


The other non-zero Christoffel symbols necessary to calculate the Riemann tensor 
components are = Tgg sin^ 6, T‘^g = cot 9, and T® ^ = —T‘^g sin^ 9. At the plane 
9 = TT12 one has 


pZ _ pZ _ 

^ ee — ^ ~ 


(6^ -I- P) (^Pe '^ + — b{P + e)y/W^n? 


el 


(32) 


and r^g = = 0, the remaining Christoffel symbols being zero or given by 

symmetry of the ones just calculated. The Riemann tensor components appearing 

(33) 


in Eqs. (27l-(29l, assuming 0 = 7r/2, are 




TdO _ _ 

^iie — ^iiip — 


2/3e i.2 ^/{b'^+l^p-P{P+e) 

2ehVP+R 


2/3Y^(b2+i2)3-b3(^+e)eS^ 




2b2(b2-|-Z2)2 ^b(/3-|-e)e -%/b2+R^e ^ 

_ 

pe y/b^+P-b(0+e) 
ey/b^+P 


(34) 

(35) 


These equation complete the geodesic deviation equations (27l-(29l. 


3.4.2. Tidal effects for pure radial motion 


Now, looking at Eqs. (27l-(29l, in general, it is not possible to find analytical so¬ 


lutions since the fiducial test particle’s velocity must also be taken into account. 
So, to solve this system of differential equations, we must first solve the equations 


of motion, Eqs. (22l-(24l, and in general there is no analytical solution. We can 


however work out some features of the tidal effects in the case of pure radial motion 
of two test particles such that throughout the motion the line that connects the two 
particles is purely radial, i.e., the case where n®, , and are zero throughout 

the particles’ motion. In this case Eq. ( [27| simplifies to 

^ = -n‘{dig^ + T\ig^) , (36) 


and Eqs. (281 and (291 yield = 0 and = 0, respectively. Although it is 


still not possible to solve analytically Eq. (361 we can instead make a qualitative 


analysis to infer the general behavior of two initially close test particles describing 
pure radial motion. 


Defining the right hand side of Eq. (36) as 


A{l) = -dig^ -T^g^ (37) 

we present in Fig. its behavior. We see that the sign of the function A changes at 
some value of the coordinate I, Ic, say. Now, when I > l^. the function A (1) is positive. 
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assume / > 0). 


(here 


A (1) > 0, which implies that the separation vector has positive acceleration. So, 
suppose two test particles describing pure radial motion, such that the line that 
connects them is purely radial. One particle has radial coordinate h, the other I 2 , 
and let us assume, h < I 2 and li ^2 > Ic with Ic > 0 for simplicity. In this case, since 
A (^) > 0, the particle with coordinate li is accelerating more than the particle with 
coordinate I 2 , and they fly apart. On the other hand, for I < Ic the function A (1) 
is negative, A (1) < 0, so the separation vector has negative acceleration. Assuming 
now li ^2 < 4 1 and still li < I 2 , this means that the particle with coordinate h is 
accelerating less, but still towards the wormhole’s throat, than the particle with 
coordinate l 2 - Numerically solving Eq. (361 for two initially close test particles in 
the regions I > Ic and I < Ic, indeed, we verify this conclusions. Notice that as 
we consider particles further away from the wormhole’s throat we recover the tidal 
behavior expected from Newtonian gravitation since the space is asymptotically 
flat. 


4. Conclusions 

A static, spherically symmetric wormhole in an enhanced Newtonian theory of grav¬ 
itation was constructed. The Newtonian wormhole’s mass density is positive, the 
gravitational held of the wormhole is well-behaved and the matter that sustains it 
has positive pressure. Test particles’ motion in the wormhole gravitational held and 
tidal effects were studied. We have also argued about the possibility of having true 
Newtonian black holes, i.e., Newtonian objects that have regions from which any 
particle must have infinite acceleration and thus infinite velocity to escape to the 
outside of it, in this enhanced Newtonian theory of gravitation. 
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